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fin We use gauge/gravity duality to study simultaneously the mass spectrum and the thermo- 

i— ' dynamics of a generic quasi-conformal gauge theory, specified by its beta function. The beta 

I function of a quasi-conformal theory almost vanishes, and the coupling is almost constant 

^ between two widely separated energy scales. Depending on whether the gravity dual has 

^^ a black hole or not, the mass spectrum is either a spectrum of quasinormal oscillations or 

*yv a normal T = mass spectrum. The mass spectrum is quantitatively correlated with the 

cn thermal properties of the system. As the theory approaches conformality, the masses have to 

f — . vanish. We show that in this limit, the masses calculated via gauge/gravity duality satisfy 

O expected scaling properties. 
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1 Introduction 

In this article, we shall consider gravity dilaton -systems dual to quasi-confornial gauge 
theories. Their characteristic feature is the presence of two widely different energy scales, 
Aetc cind AtC) between which the coupling constant of the theory evolves very slowly; the 
theory is almost conformal, the beta function almost vanishes. Below Atc a-nd above Aetc 
the coupling runs similarly to asymptotically free theories. 

Quasi-conformal gauge theories are of current interest for beyond the Standard Model 
physics in the context of walking technicolor [Ij. Our model lacks proper flavor dynamics 
of technicolor, but this connection nevertheless motivates our notation and the quantitative 
estimate Aetc ~ 10'^ Atc- 



Our gravity dual, Eq. (2.2) below, is based on the improved holographic QCD presented 



in [21 [21 S] and further developed in O El [7|. The extension of the model from QCD to 
technicolor is based on building the required physics in the ansatz for the beta function of 
the coupling g^ of the theory, in a certain scheme. For this we use the ansatz 
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:i-A)^ + e 
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which near A = 0, 1 and oo can be expanded or approximated as 
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In these one has built asymptotic freedom in the ultraviolet (UV) domain, i.e., as A — t- 0, 
walking near A = 1 provided e is small and confinement in the infrared (IR) large A domain 
if a = 2c/3 [2J. The values of parameters c,a,e reflect the physics one wants to describe: 
c>10 fixes the high T phase transition to be of first order [6] and e has to be a small number 
to describe near conformality. The case e = leads to an infrared fixed point (IRFP) [7J. 

Thermodynamics of the model has already been studied in [BJ: the system can exist in 
three different phases with phase transitions at T = Ttc ~ Atc and T = Tetc ~ Aetc • The 
main goal of this article is to study the spectrum of the theory, which is either a usual T = 
mass spectrum, if there is no black hole in the gravity dual or a quasinormal spectrum ^ with 
complex energies, if there is one. The relation between the numerical values of TtC; 2etc 
and the mass values is then also obtained. 

Of particular interest will be the behavior of the mass spectrum when conformality is 
approached, e — )• 0. We shall show that in this limit the mass spectrum as computed from 
gauge/gravity duality satisfies Miransky scaling [5] 



m,(c, e) = A{c^ e) exp 



2 1\ TT 

3 c) y/e 



:i.5) 



where A{c,e) is slowly varying when e — t- 0. 



The bottom-up gravity dual presented here lacks detailed flavor dynamics and is thus 
not intended to be holographic description of walking technicolor. Examples of top down 
technicolor duals are given in |101 111) . In any case its merit is that it presents a concrete 
gravitational dual with phase transitions at two widely separated energy scales. A general 
discussion of these in terms of multiple big black holes is given in Section G of |3j . 

Section 2 discusses the gravity dual and, in particular, the determination of the dilaton 
potential which, after solving the Einstein equations, leads to the beta function ( 1.1 ). Details 



of numerical evaluation of the background are also presented. Section 3 reviews the thermo- 
dynamics, already discussed in [6]. Depending on parameters, the model can describe either 
first order or continuous transitions. General principles of determining the masses and quasi- 
normal modes from Green's functions are summarised in Section 4. One practical limitation 
of the present dual is that it can only be computed numerically. However, one can present 
a simplified version of the model containing qualitatively its essential properties [20], this is 
applied to the determination of the quasinormal modes in Section 5. The full computation 
of the spectrum is carried out in Section 6. 



2 Gravity dual 

2.1 Equations 

The gravity equations of the model are as follows |3l|6]. The gravity action, including the 
scalar field (j){z), in standard notation and in the Einstein frame, is 



s = Ye^J d'-V^9[R- lid.^f + v 



With the metric ansatz 



ds^ = b'^{z) 



-f{z)dt^ + dx2 + 



dz 



2 n 



/(^) 



(2.1) 



(2.2) 



the three functions 6(2), f{z) in the metric and the scalar field 0(z) are determined from the 
three equations (5 = 6'(z), etc.) 



62 h bfb'^,, 

f b 
4+3- = 0. 
/ b 

Further, from the functions so evaluated, the beta function follows as 



(2.3) 
(2.4) 
(2.5) 



/3(A) ^5^ ^5^, A(.)^e^ 



~ g^N,. 



(2.6) 



Thus A = e^^^> is the coupHng at the energy scale h{z). Often it is practical to use the rescaled 
form 

«^' (2^7) 



x = 



3A 



and to define 



W 



b ■ b"^ b 

1?' "'' = ■'»- 1- 



(2i 



which changes (2.4) to first order form. Inserting further W = W'{\)\ and then A = 6/3/6 



-bWP permits one to solve W{X) in terms of /? (Eq. (2.11) below). 
Note the two roles, mathematical and physical, played by /3(A). Irrespective of the physical 



context this combination of A, 6, arises very naturally in the numerical solution of (2.3)-(2.5) 



(see, e.g., (2.9) or (2.10) below). Physically it models the running of the coupling A(z) as a 
function of the energy scale ~ b{z). 



Note further that /3(A) can be solved from (2.6) only after all the fields have been solved. 

le cai 



This requires a knowledge ofV{X). However, given /3(A) one can from (2.3) with / = 1 solve 

2' 



Vo{X) = 12W\X) 
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3A 



1. 



(2.9) 



Here W{X) can be solved in terms of the beta function (2.6) and, introducing the explicit 
ansatz (|1.1|), one gets 



W{X) 
VF(0) 



exp 



\-x^^^^ 



(2.10) 



exp< WE 



2\/3a^/^(-2 + a^/^(l + e))(arctan[ '^^ ] - arctan —=)+ 



ai/3(2 + ai/3(i + e)) log[(l + a^/'^Xf/{l - a^'^X + a^/^A^)] + 2 log[l + aA^])] } , (2.11) 
where VF(0) is the inverse of the AdS radius: 

P^(0) = ^. (2.12) 



1 



At large A, Eq. (2.11) implies that W{X) ~ A^^/^^") so that Vq{X) ~ A^^/^^'^), more detailed 



asymptotic expressions are given in [20] and below in Eqs. (2.23) - (2.28). As shown in [2], 
confinement requires that the large A limit he V ^ A^'^(log A)*^^''" with a > 1. Thus we 
fix 

a='^c. (2.13) 



A mass spectrum of type rn? ~ n ^2j further implies a = 2. This is also built in Eq. (2.27) 
below. To avoid spoiling the UV small-z, small-A behavior we shall use the potential 



V{X) = yo(A)i 



4og(F + A^) 
logF 



(2.14) 



Here F is a parameter which sets the scale at which confinement effects set in. Since the 
quasi-conformal effects operate, by construction, at A ~ 1 we certainly expect F^'^ ^ 1 and 
choose F = 1000, to have a large separation between ETC and TC scales. For comparison, 
in |3] the potential for SU(A'^c) Yang-Mills thermodynamics was chosen to be 

^^ihqcd(A) = ;§{i + 0.04128A + (7.35A)^/Vlog[l + (13A)2]|. (2.15) 



We re-emphasize that V{X) in (2.14) is given in terms of a quantity /3(A) which one knows 
only after one has solved the problem. So one has to numerically check that the output beta 
function to sufficient accuracy reproduces the input one; this seems to work very well [B]. 

2.2 Numerical solution 



The numerical solution of Eqs. (2.3) - (2.5) with V given by (2.14), means finding a set of 
functions b{z) , f (z) , X{z) parametrised by the horizon value A/j = X{zh), f{zh) = [H E] • 
Various physical quantities computed from the set of functions are then functions of A/j. For 
example, see Fig, [g] for T = T(A/i). 




To solve Eqs. (2.3) - (2.5) it is practical to use the W in (2.8) to write them in the form 

^{Wf + lbV), (2.16) 

(2.17) 
(2.18) 
/ = 3fbW, (2.19) 

Details of numerical integration can be found in Sect. 2.2 of [6]. We add here a number of 
points relevant for the purposes of this article. 

1. The choice of the set of Xh values depends on whether one studies thermodynamics 
or the mass spectrum. In the former case one has to integrate over a large range of 
temperatures so that one chooses a large set of A/j values from small values, < 1 for the 
UV and ^ 1 for the IR. For the mass spectrum one should find a special no-black hole 
solution f(z) = 1. How this is to be done is discussed in Appendix A. However, one 
can also possible to start from a / 7^ 1 numerical solution with one very large value of 
Xh and correspondingly very large value of Zh- Staying within the values of < z < Zh, 
for which f{z) is very close to unity, one has an approximate solution of b{z), X{z) for 

/ = !• 

2. The initial condition bh = b{zh) follows from b(zh) = exp[J '" dX/ (3(X)]: 

cil + e)logbh = —--^log[c{l + e)Xh] 
Xh 1 + e 

l-e + a(l + e)2/ 1 - A;, ^ l\ 

-\ i=r- r arctan ^ arctan — ^ 

Ve(l + e) V Ve VeJ 

2 1 + e) ^ 1 + e ' ^ ^ 



where the constant of integration is fixed so that the leading UV A — )• behavior is the 



2-loop one, i.e., the two last terms in (2.20) vanish at A/i = 



3. The initial condition Wh = bhV{Xh)/{—'ifh) follows from (2.3) by requiring that all 
fields remain finite at the horizon. For f^ one can take an arbitrary negative value. 
The energy unit is determined so that at z — t- 

>/ N 1 ^1 loglogil/fAz)! , , X , . s 

where we choose units A = 1. 

4. Due to large ETC and TC scale separation large values of z enter into the analysis and 
it is useful to have an analytic approximation to numerics. One cannot analytically 
solve the equations at large z and large A, but one finds that a good approximation for 
b{z) is[] 

b{z) = boe-^-''>\ (2.22) 

H here is an IR scale, analogous to the A in the UV. Another definition will be 
given below in ( |6.2[ ). From the numerical solutions we can determine S, and we find 
that typically H/A ~ 0.001. This is consistent with the large energy scale splitting 
ETC/TC~ 1000. Inserting this to the equations of motion above one finds that, up to 
corrections 1 + 0{l/z^) ~ 1 + 0(1/ log A), 

b{z) = boe-^~'^\ 
W{z) = |e(=^)'(Hz), £ = ^, 

m = -iA(l + ^ + C)(^) 
V{z) = 9W^iz) = §e^^-'y\Ezf = 

z 
-f{zh) = A^T = Qr?Zh. (2.28) 

Here A = A/Aq, where Aq is the constant of integration. These equations will in practice 
(Appendix A) be important as the initial conditions for computing the / = 1 gravity 
solutions. 

5. As a concrete example of computed bulk fields, Fig. [T] shows zh{z), \{z), W{z) for 
one set of parameters; Xh = 10^ with Zh = 1805 was used here. This also shows 
how the required structure with two widely separated energy scales is refiected in the 
bulk fields. For example, the coupling \{z) first approaches an almost stable IR fixed 





(2.23) 




(2.24) 


xf'\ 


(2.25) 
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(2.26) 


§A^/3(ilogA)V2, 


(2.27) 



^ Note that one may use zh{z) = 6oexp[— (H^)^] to capture both the UV and IR behavior, as in the model 
of Section p] |20]. For this discussion of the behavior deep in the IR, it suffices to take only the Gaussian tail 
into account. Formulas for a more general parametrisation 6 ~ exp[— (Hz)"]2;^ are given in [20j . 
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Figure 1 : The functions zb{z), A(z), W{z) for e ^ 0.1, c = 13/(1 + e), F = 1000. Tlie daslied curve 
is the fit (2.29) with parameters given in the text. The liorizon is at Zh ~ 1805 and at this point A(z) 



has grown to the value 10 



point A = 1 near z = 1. However, instead of freezing at this value, A gradually starts 



growing and ultimately at large z grows as specified by equation (2.25) valid at large 
z: A ~ exp(| H^z^), W ~ exp(S^2;^). The transition between a UV fixed point and a 
quasi stable IR fixed point is particularly well seen in the behavior of b{z), which first 
behaves as Cu/z, then as Cj/z with a decreasing AdS radius and ultimately in the 



far IR shows the gaussian decrease (2.23). In fact, the numerically computed zb{z) in 



Fig. [T] can be fitted by an expression of the type 

Z b{z) = [Cu - {Cu - Cl)Q{z - Zqir)] e"^'^' , 

where in Fig. [l|£[/ = 1.1, £/ = 0.18, G(z-Zqir) = \+\ tanh((z-0.6)/0.8), H = 1/730. 



(2.29) 



3 Thermodynamics 

The thermodynamics of the present model has been discussed in [6j, but with a different 
functional dependence on the parameter F in the dilaton potential. More precisely, in [B] the 
F-dependent term in (2.14) was chosen to approach one as e — )• 0, while now the confinement 
factor in (2.14) is independent of the walking parameter e. This is more realistic since 



confinement enters at large A, far from the IR fixed point A = 1. 
To summarize the results, we shall choose as reference values 

e = 0.1, c(l + e) = 13, F = 1000. 



(3.1) 



The combination c(l + e) is the leading Hoop coefficient of the beta function ( 1.1 ). One then 
generates numerically a family of solutions 6(z), X{z), f{z), f{zh) = 0, parametrised by the 
value Xfi = X{zh), and normalises the solutions as in (2.21). From the solutions one then 
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Figure 2: The function r(A/,.) obtained for e == 0.1, c = 13/(1 + e), F = 1000. The first minimum is 
at Xh ~ 0.218, T/A = 0.8 the second at A^ — 20, T/A = 0.001. The corresponding transitions are at 
T/A = 0.653 (Fig.JS]) and at T/A = 0.00098 (Fig.|4|- The behavior at large A^ (inset) can be fitted 
by 0.00050 log A^■'*^ 
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Figure 3: The equation of state obtained for e = 0.1, c = 13/(1 + e), F = 1000 for the region 
around the first order ETC phase transition at T^tc = 0.602A. Dotted lines are unstable, dashed 
lines metastable (supercooled or -heated) . The boxed confinement transition to the low T phase with 
p = at Ttc = 0.000977A is blown up in Fig. ^ 



computes thermodynamics by evaluating 

AttT = -f{zh), 

I rco _ 

p = -v-^ I dXh 



dT 



4G5 Jx, 
e = Ts — p, 
^2 _ s{T) _ /3(A) dT 
^' ~ Ts'{T) ~ ~3T~ dXh 
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Figure 4: The equation of state obtained for e = 0.1, c = 13/(1 + e), F = 1000 in the region around 
the confinement transition to a low T phase with p = (thick red line) at T-pc = 0.000977A. Dotted 
lines are unstable, dashed lines metastable (supercooled or -heated). 



Each quantity can be regarded as a function of Zf^ or X^ or T. If the extreme UV pressure is 
5cffvr2r^/90, then 

Thermal physics is basically contained in the function T(Aft), shown in Fig. pi The 
model is constructed so as to have two transitions and, accordingly, there are two minima. 
The mechanically stable regions with c^ > correspond to dT/dXh < (cf. Eq. (3.5) and 
remember that /3 < 0) . With the structure of Fig. ^ one can have T equal in two different 
decreasing branches of T(Xh), and if they also have the same pressure, one has two different 
phases in thermal and mechanical equilibrium, i.e. a first order phase transition, see Fig. [3| 
The transition temperature Tetc is somewhat larger than T at the first minimum of Fig. [2| 
In view of the resemblance to walking technicolor, we call this the extended technicolor 
transition. 

If one continues towards lower temperatures, at some T the pressure must go to zero, see 
Fig. |4J It is zero at Xh = oo, as defined by Eq. (3.3), and since dT/dXh > 0, it is negative 



over some range in the IR. There is a cusp in p/T^ at the larger minimum of T{Xfi) in Fig. pi 
Below this temperature, p starts growing and passes zero at some T = Txc- The other stable 
phase is the "hadron gas" phase with p = 0, in this strongly coupled, large Nc theory, and 
the transition is the analogue of the deconfining QCD transition. The large scale separation, 
^ETc/^TC ~ 600 is due to the choice of a large value for the parameter F in (2.14). 

A feeling of the numerical values involved, and on the general sha pe of T{Xh), can be 
obtained by adding the UV behavior f{zh) ~ ^/zh and the IR limit in (2.28) resulting in 



ttT 



Zh 



+ 



Zh 



Aei/(^"^'^)(6oA,)^^/''o + 3H(|logA,)i/2. 



(3.7) 



Thus there is a very rapid decrease from the UV and a very slow increase in the IR. This 
is also reflected in the numerical value of the IR parameter H/A ~ 0.0014 in (2.29). The 
fit 0.00050(log A/i)*^'^^ to the curve in Fig. ^is compatible with this value of H. We remind 
that the power 1/2 of log A/^ was imposed by choosing the confinement factor ~ (logA)^'^ in 



(2.14). 



4 Spectrum of the theory: masses and quasinormal modes 

The mass spectrum is obtained as poles of a Green's function of an operator with appropriate 
quantum numbers in the background (2.2) with f{z) = 1. The poles corresponding to stable 
particles are on the real axis; if the backround has f{z) with a horizon, f{zh) = 0, the poles 
move to the lower complex u plane. 

To determine this Green's function in gauge/gravity duality we need the equation satisfied 
by fluctuations in the background (2.2). The fluctuations h^iy of the metric satisfying d^h^^ = 



h^ = 0, i.e. the tensor fluctuations, are simple since they do not mix with other fluctuations. 
For dilaton fluctuations, i.e. the scalar fluctuations, one must take into account that they mix 
with the fluctuations of the metric. The appropriate equation has been derived in |13| 114] . 
Keeping flrst / 7^ 1 and taking fe = the equation for the gauge invariant scalar fluctuation 

(h = Moj, z) is 



+ ^log(/53).f 



X 
X 



d 

dz 



2u3\ 



logif'b 



X 
X 



+ 



UJ 



0, 



(4.1) 



where X was deflned in (2.7). For the tensor fluctuation one simply removes the X-terms. 
For the true T = mass spectrum one takes / = 1 and replaces uj"^ 



UJ 



k' 



m 



m- 



The methods of computing the Green's functions and their poles from (4.1) are standard 
|16 1ll7l[T8l ll9j. To interpret the outcome it is convenient to transform (4.1 ) to a Schrodinger- 
like form with cj^ as the eigenvalue. In order to achieve this, one removes the 1//^ in (4.1) 
by introducing a new variable 

dz 



u 



which transforms the interval 0<2;<z/jtoO<u<cxD. Introducing further 



the equation becomes 
where 



^p{u) = vh^ 4){u 



r(S) 



i;"{u) + Vf\u- ZhMu) = u;V(w) 



V, 



(S) 



T'{u;z^) = P 



36 352 3/6 X /36 2f\X 
26 462 2/6 X V ^ / y ^. 



(4.2) 

(4.3) 
(4.4) 
(4.5) 



z=z{u) 



The superscript S indicates that we are considering scalar fluctuation s. F or tensor fluctuation 



r{T) 



the analysis is similar, and the result for V) (n; z^j is obtained from ( 4.5 ) by simply removing 



the X-terms. When u w z — )• the requirement b{z) — )• C/z impUes that the potentials behave 
as 



and when u — )• oo, z — )• Zh, 



where 



4^2 

rii) 



Vf\u;zh)^A('\zh)ef'^-^0, (4.7) 



When f{z) = 1, the Schrodinger-type equation becomes 



'Sch 

where 



r{z) + V^i{z)i;{z) = m'i;{z), i = S,T (4.9) 



36 36^ X bX 



AS) 

26 ' 4^2 ' X ' "6X 



^Sch^ = 7^ + 7I^ + ^ + 3i^. (4.10) 



(T) 

For the tensor Schrodinger potential, Vg^h i the X terms are omitted. For a simple estimate of 
the IR behavior, one may take for b{z) the numerically accurate fit (2.29), zb{z) = exp(— H^z^^ 



valid for large z. This leads to the approximate confining tensor potential 

^sT = ^+6-^ + 9-'-^- (4-11) 



This has the minimum value (3vT5 + 6)5^ at H^z^ = y^5/12 and the mass spectrum m^ 



12H (n+2). However, the existence of the dynamical ETC and TC scales produces interesting 
structures in between. 

5 Spectrum in a model of the deconfinement transition 



To appreciate the relation between / = 1 and / / 1 cases in the dilaton background (2.2), 
we shall first consider a simplified version of the TC or deconfinement transition [2Qj . It is 
based on the ansatqj 

biz) = -^e-h\ (5.1) 

from which one can derive A(z),/3(A) and, in particular. 



z 



l)e< + 1 '^^ zl-l + e~'h 



^- x;c. .. -r X -■■ ^^ 



The tensor potential in this case is 

0^)^1 + 1-^+2..- (5.3) 



^For simplicity, we fix the IR scale 'E? = 1/3 in this section. 
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given by (4.51 with X-tcmis omitted for various 



Figure 5: Left panel: The tensor potential V^ 

(T) ' ' (T) J I 

values of Zh- For increasing Zh the potential V^ approaches the potential Vg(,j^ in (5.3). Right panel: 
The corresponding quasinormal tensor frequences computed from (4.1). The mass spectrum (5.4) is 
plotted on the real axis. Poles with Re (w) — )■ — Rew are not shown. 



and from this the tensor glueball mass spectrum becomes [12j 

2V2 + n, n = 0,1,2,.... 



rrir, 



The curve T = T{zh) has a minimum Tmin = 0.3862 at Zh 
a high T deconfined phase for T > Tc = 0.4000, z^ < Zc 



(5.4) 

1.466 and the system is in 
: 1.30. Inserting (5.1) and 



(5.2) to (4.10) and performing the change of variables (4.2) numerically, one finds the tensor 
potential plotted in Fig. ^ The overall pattern is obvious: For small n ~ z, f{z) = 1 and 
Vf {u; Zh) ~ V^ll^{z), both for scalar (i = S) and tensor {i = T). The effects of Zh start being 
felt and f{z) starts to decrease towards f{zfi) = when u ^ z approaches z^. Furthermore, 
as z approaches z^, the variable u approaches infinity and, due to the f'^{z) factor in (4.5), 



the potential V^ 



W, 



/ 



CO] Zh) approaches zero. 



The mass spectrum (5.4) can be derived by solving the Schrodinger equation (4.9). To 



solve the quasinormal spectrum when / 7^ 1 it is actually most convenient to go back to (4.1 ) 



and use standard techniques for computing the Green's function. To this end one solves first 
Eq. (4.1 ) numerically by starting the integration at z 



Zh expanding for small z — ZhSO that 
Secondly, one constructs analytically the two linearly 



the leading term is ~ (z 

independent solutions, 0„ ~ 1 and (/>„ ~ 2;* near z = by expanding the solution at small z. 

Finally, the quasinormal modes are those values of w = Re u + zim uj for which the solution 
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e = 0.1,c = 13/(l+e) 
F = 1000 




100 1000 ^ 10^ 



Figure 6: The scalar (dashed hne) and tensor (continuous line) potentials in the IR large-z TC 
region for e = 0.1, c = 13/(1 + e), F = 1000 (note logio scale). The large-z tail is given by (4.11 1 with 
^ sa 0.0014A. Two lowest scalar and tensor excitations are p lotte d, the ordering is £'q < -Eq < 



AS) 



AT 



E\ < El ' . The corresponding mass values, m, = E, are in (6.11 



~ z^ at small z and the coefficient of (f)u vanishes. Examples of the lowest quasinormal states, 
plotted in Fig. [Sj are 






4.0-1.78i, 6.36 -3.3i, 8.74 - 4.78i 

2.94-0.32i, 4.13-0.77i, 5.46 - 1.22i, 6.8 - 1.66i 

2.83 -O.z, 3.46-0.004i, 4.02 - 0.036i, 4.57-O.lli, 5.19 



0.20i 5.85-0.29i 



The pattern is clear: for large Zh the quasinormal frequencies approach the real mass spectrum 



(5.4). However, then the system is in the unstable small black hole phase with T'{zh) > 0. 
With decreasing Zh the spectrum moves away from the real axis, and approaches uj = 2n(ibl — 
i)/zh = 27rTn(ibl — i) as the system enters the stable high T phase. For example, Zh = 1 
corresponds to T = 1.08Tc, but then the states are very broad as they are melting away. 



6 Spectrum for quasi-conformal theory 

To have an overall picture of the situation, we first plot in Figs. m\ (on logar ithmic scales) 
and y (linear scales, large V^^^) the / = 1 Schrodinger potentials in (4.10) choosing the 
same reference values (3.1) as in the discussion of thermodynamics: e = 0.1, c = 13/(1 + e), 
F = 1000. Later we shall study the approach to an IR fixed point by taking e — )■ 0. Note 
that formally the potentials are plotted in units of A, Vgcj^/A^ vs Az, but we have chosen 
A = l. 
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Figure 7: The scalar potential and the energy {Eo — 0.355 
scalar state (note hnear scale). 



IQ- 



and wave function of the lowest 



Overall, the potentials follow very well the pattern in (4.11) with H ~ 1/730 as in Fig. [Tl 



with some important modifications: there appear local minima at the ETC scale and the 
scalar potential has more structure than the tensor one and also binds more strongly. 



6.1 Stable states and quasinormal modes at the TC level 

The lowest energy stable states of the system are those near the minimum of the potential 
^Sch- Computing the four lowest scalar and tensor stable states one finds that their mass 
values are (as always, in units of A = 1) 



m 



m 



(S) 
(T) 



(0.005948, 0.0077235, 0.0090700, 0.010182) 
(0.0069018, 0.0083675, 0.0095586, 0.010581) 



(6.1) 



In spite of the prominent dip in the scalar potential the lowest scalar mass, the analogue of 
standard model Higgs, is thus not particularly well separated from the higher excitations. 
Its wave function is plotted in Fig. [7| and one sees that the narrowness of the potential well 
pushes the state higher. Masses should be compared with Ttc = 1-0 x 10~^A, see Fig. [4| and 
one observes that the lowest mass ~ 2ttT'j-c. 



The tensor potential is well approximated by the form (4.11). In fact, one could use this 
to define the IR scale H via the equation 



:(c,e) = ^^/i^, V'{Zrain) = 0- 



(6.2) 



Numerically, this is close to the alternative definition (2.23), b{z) 
0.1, c(l + e) = 13 the former gives E = 1/680, the latter E = 1/730. 



exp(— H z ); for e 



Since the approximation (4.11) works so well, the approximate quasinormal tensor modes 
can be read from the results of subsection [sj by noting that there 3H^ = 1, i.e., in Fig. [sj 
one should replace to — )• \/^E. For the scalar case the quasinormal modes have not been 
computed. 
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6.2 Stable states and quasinormal modes at the ETC level 

The / = 1 potential over the entire range was plotted in Fig. [6| A close-up of the ETC 
region, 2; ~ 1, A ~ 1, V5 ~ 1, is shown in Fig.jSj left panel. Overall, the potential is confining 
but extremely wide, one has to go to 2; ~ 10^ before the large z part reaches the magnitude 
100 of the peak in the ETC region. Due to the confining nature of the potential, there are 



stable states of large excitation number. In fact for the potential (4.11), which approximates 
very well the potential in Fig. [6j one has 

m^ = 12H2(n + 2). (6.3) 

Thus, for m = 1 the excitation number would be n = 44410 and the mass splitting Am = 
y^3/nH. However, there also are local minima in the potential. These arise due to the 
quasiconformal behavior of the beta function near A = 1 and the related rapid variation of 
the bulk fields in Fig. [Tj 

Potentially, the most interesting feature is that the scalar potential even dips clearly below 
zero, so that even stable states are thinkable. However, for all the parameters the dip seems 
to be so narrow that the energy levels are lifted to the metastability region. 

One may ask what the basic reason for the appearance of this structure in the ETC region 
is. It really comes from the transition of 6(z) from the UV conformal region with 6(2:) = Cu / z 



to a quasi-conformal region with b{z) ~ Ci/z. This is described by Eq. (2.29) and inserting 
this to Vsch one sees that the peak structure is formed. 

To see the effect of minima at F > 0, forget first the increase of the potential at very large 
z. The minima then would correspond to metastable states. To estimate their energies, one 
may demand that the wave function of the potential in Fig. |8]vanish at z = 0.01 and z = 1.2. 
One finds that this potential can bind two states with masses 

mo = 4.42, mi = 8.92, (6.4) 

shown as horizontal lines in Fig.|8j Of course, these states are only metastable. The tunneling 
actions for these states are of the order of 1 so that estimates of widths are are very unreliable. 
Restoring the increase of the potential at large z one sees that these states simply disappear 



in the dense spectrum of large m states (6.3), effectively doubling two of the states (6.3). 



However, the levels (6.4) will reenter when considering quasinormal modes. 



The appropriate / / 1 potentials are plotted in Figs. [8] and [9j Fig. [8] compares the / = 1 

/ 



potentials V^ll^{z) and / / 1 potentials V^^ {u; Zh) at T = Ttc- Fig. 9 left panel, shows 



the / dependent potential V^ {u;zh) for a wider range of z^. To appreciate the relation 
Xh = X{zh), see figure [Ij One sees that, as suggested by Fig. [sj when z^ approaches the UV, 
the metastability peak is washed away. 

The quasinormal spectrum computed for various values of Zh is also shown in Fig. [9} For 
large values of Zh (A/t>l the metastability peak is efficient and there are states with negligible 
imaginary part. These are the metastable states suggested physically by the existence of the 
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Figure 8: Left panel: The scalar and tensor potentials in the UV small-z ETC region for e = 0.1, 
c = 13/(1 + e), F — 1000 (note linear scale). For z — >■ both potentials approach ~ 15/ (4z^). For 
z — > oo one has to go to z « 10^ before VgV) grows back to the value 100 (Fig. 6|. Right panel: 
The / 7^ 1 scalar and tensor potentials Vf {u;zh) in (4.51 for the same values of e,c,F and for 
Xh^Xc^ 0.71, zh = 0.71, T = Ttc = 0.6A. 
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A,, = U.Tn 
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-Ar^Jv329- 






( s) 
Figure 9: Left panel: The dependence of the scalar potential V} {u,Zh) for A^ = 

0.138, 0.329, 0.520, 0.711, 1.0, Zh = 0.57,0.719,0.724,0.92,2.96, ttT « l/z^. The peak disappears 

when one is approaching the UV, when Xh decreases. Right: Corresponding quasinormal modes. 

peak in the potential. Concretely, the lowest states for four of the X^ values are 



Xh = l,zh = 2.96 
Xh = 0.711, Zh = 0.922 
Xh = 0.520, Zh = 0.724 
Xh = 0.329, Zh = 0.719 



4.41 - 0.0036i, 8.89-O.OlOi, 12.9 - 0.07H, 17.9 - 0.069i 
4.33-0.018i, 8.86-0.18i, 13.5 - 0.847i, 21.4-2.50i 
3.44-0.26i, 7.90-0.79i, 12.0-1.68i, 20.2 - 3.62i, 
0.00-1.17i, 6.32-2.09i, 10.2 - 3.47i, 29.5 - 4.80i. (6.5) 



Comparing with (6.4) one sees that for A/i>l the temperature is already so low that the 



real parts are very close to the estimated energies of metastable states. However, now one 
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also has a controlled determination of the imaginary parts. For increasing vrT ~ l/-z/i the 
imaginary part rapidly grows, the states are melted away. The real part of the lowest state 
starts decreasing and soon is even very close to zero. Note that this is quite different from 
what was observed with the model computation plotted in Fig. [5] The reason is obvious: 
the pattern in Fig. [9] is determined by the effective vanishing of the ETC metastability peak 
when T grows, while there is no such effect in the TC region. 

Finally, we consider how our results are affected by the choice of the parameters in the beta 
function, i.e. in the dilaton potential. First, since for e 7^ the theory always confines at 
large A, the large z behavior and the spectra at the TC level are practically independent on 
the choice of parameters as long as e 7^ 0. Second, at small z, i.e. at the ETC level, there arise 
structures more sensitively dependent on the values of the parameters. This is illustrated in 



Fig. 10 which gives examples of how the Schrodinger potentials vary when parameters are 



changed. In the plots, we show the potentials scaled with the IR scale S defined in (6.2) vs. 



Sz. As we have already noted, the shape of the potential at small z consists of a minimum 
followed by a peak as z increases. There are two main effects which arise as the parameters 
in the potential are changed: The first effect is that the location of the high temperature 
minimum shifts towards larger values of z as the parameter e is increased. The left panel 
of the figure shows how this happens for two values e = 0.1 and 0.3. The parameter F 
has similar effect, i.e. larger value of F implies larger separation between the high and low 
temperature minima. The second quantitative feature is the depth of the high temperature 
minimum and the height of the peak. This is mainly controlled by the parameter c as can be 
seen from the left panel of Fig. 10 and in more detail for e = 0.01 from the right panel. In 
conclusion, the qualitative shape of the Schrodinger potentials and our corresponding results 
for the spectra appear robust as the parameters are varied so that the walking nature of the 
theory is preserved. Of course the quantitative results like the number of metastable states 
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Figure 10: Left panel: The Schrodinger potentials scaled with the IR scale (6.2) for two values of 
e, plotted vs Sz. Right panel: Dependence of the Schrodinger potentials on c at fixed e = 0.01. If 
one increases e to 0.1, the height of the peak and the depth of the dip at c(l + e) = 14.3 increase by 
~ 30%. 
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at the ETC level, and how rapidly they melt away at increasing temperature will depend on 
the exact numerical values chosen for the parameters. 



6.3 Approaching the IR fixed point, e — )■ 

When e — )• in (2.2), the /3 function approaches one with an IR fixed point at A 
remind that, in practice, always a = o c): 



/3(A) 



1 + a 



[{l-Xf + e]. 



1 (we 



(6.6) 



In the limit the theory is conformal both in the UV and in the IR. Thus one expects the 
mass spectrum to approach zero, in the potential language the confining part +9H'^z^ of the 
potential disappears. How this happens can be derived as follows. 



Integrating dX/dt = /3(A), t = log(/i//io) in the approximation (6.6) one has 

A(t) = 1 - Vetan(— ^Vet), 
1 + a 



(6.7) 



normalised so that A(0) = 1. This decreases monotonically from +oo, reached at the value 
i < defining an IR scale: 

/iiR (l + a)7r 



tiR = log o /- • 

The UV scale can be defined as the value of t > for which X{t) has decreased to zero: 

(l + a)7r 1 + a 



*uv 



2cV^ 



-tm- 



(6.8) 



(6.9) 



Thus 






e*iR-*uv = exp 



2 1 
3+c 



vr 



(6.10) 



where a = | c was inserted. We thus expect that, at small e, the IR scales and masses behave 
as 



m 



^exp 



D 



D 



2 r 

3 c 



(6.11) 



A = constant. 

This prediction works surprisi ngly well. All the masses are proportional to the IR scale 
H = (5/12)-'^/^/zmin, defined in (6.2). Here Zmin is the value at which the tensor potential 



has its IR minimum. Its dependence on e at fixed c(l + e) (keeping this fixed keeps the 



leading UV behavior unmodified by the variation of e) is shown in Fig. 11, left. For masses 
the behavior is very similar: 



m[^\c{l + e) 



m. 



(T) 



m 



iS), 



m, 



(S) 



(c(l + e) 
c(l + e) 
(c(l + e): 



= 9,e) = 5.9874 exp[-2.4348/Ve], 

= 9,e) = 6.6283 exp[-2.4337/Ve], 

= 9,e) = 7.7630 exp[-2.4342/Ve], 

17, e) = 14.671 exp[-2.2764/Ve]. 



(6.12) 
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Figure 11: Left panel: The dependence of the IR scale S defined in Ec^. 6.2 on e at fixed c(l + e) = 
9 (A = 1.4769, B = 2.4358) and c(l + e) = 17 (A = 3.60638, Z? = 2.277). All masses behave very 
similarly. Right panel: The dependence of the tensor and scalar (crosses) masses squared on c at fixed 
e = 0.1. 



According to (6.11 ) the calculated coefficient in the exponent is £> = 77r/9 = 2.4435 for c = 9 
and D = 377r/51 = 2.2792 for c = 17. 

The pre-exponential factor in (1.5) is non- universal. One finds that the c dependence 
of H or masses at fixed e is powerlike, ~ c^'*^^, with p(e) = 1.287/e'''^^* over the interval 
0.01 < e < 0.1. An example at fixed e is shown in Fig. [TTl right. 



Comparing (6.12) with (6.1) one sees that the relative hierarchy of the states remains 



unchanged when e — )• 0. One could have hoped the lowest scalar state to separate from the 
others in this limit |21] . but the model clearly needs new dynamical input to accomplish this. 



Since "Miransky scaling" (6.11) generally arises as one approaches a transition into a 



conformal theory, the prediction (6.11) and its verification by the full numerical analysis, 



provides another successful test for the consistency of the gauge/gravity framework considered 
in this paper. 

7 Conclusions 

In this paper, we have considered five dimensional dilaton gravity solutions which, in the 
framework of holography, are dual to a quasi-conformal gauge theory. In these theories 
the beta function almost vanishes, and the coupling is almost constant between two widely 
separated energy scales. 

The dilaton potential V{4') was determined so that using it as input to Einstein's equations, 
the required beta function was obtained as output. Computing numerically solutions with 
and without black hole, we determined the thermodynamics of these theories as well as 
their vacuum mass and finite temperature quasinormal spectrum. The gravity solutions 
imply that the corresponding field theory undergoes two phase transitions. Starting from low 
temperatures, there first is a transition from confining to quasi-conformal phase at Ttc and 
then, at significantly higher temperature Tetc ~ IO^Ttc a further transition to a strongly 
interacting plasma phase. 
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We determined the spectra of scalar and tensor states and studied how they are affected 
by quasiconformality between widely separated energy scales. The vacuum (no black hole 
in dual) mass spectrum is little affected, there are the usual vacuum states with masses 
m? ~ ^^TC ^P ^° very large values of n. However, when approaching conformality at 
e = 0, the overall scale of the mass spectrum was shown to exhibit Miransky type scaling, 
m ~ exp [—D/^/e], with known D. This shows the robustness of the holographic scheme: the 



approach to conformality is built in the ansatz (1.1), but the masses are computed within a 



rather complicated scheme of choosing an ansatz for the bulk metric, solving the bulk fields 
from Einstein's equations and solving field equations in this background. Still the Miransky 
scaling comes out numerically with great accuracy. Another sign of the robustness is that the 
Schrodinger potentials computed for scalar states contain negative parts which, in principle, 
could predict stable states at both Tetc and Ttc- However, it seems that these dips never 
can bind any states. 

In contrast to the vacuum spectrum, the quasinormal spectrum (black hole in the dual) in 
strongly affected. In fact, there are new states with small imaginary part at the Tetc scale, 
see (6.5). The existence of these states is already suggested by dips in the vacuum potential. 



Dynamically the origin of these dips is in the transition of the system from one conformal 
phase in the UV {b{z) = Ly jz) to a quasi- conformal phase {h(z) ^ Cj/z). 

The emphasis of this article was on quasi-conformality between two widely separated en- 
ergy scales. This is what walking technicolor is built on, but it is obvious that for realistic 
applications to technicolor models, one would need to take better into account the dynamics 
of fiavor degrees of freedom. 



Acknowledgements. JA thanks the Magnus Ehrnrooth foundation and TA the Vaisala 
foundation for financial support. 

A Calculation of the Schrodinger potentials and their eigen- 
values 



In this Appendix, we shall discuss the integration of the Einstein's equations (2.16)-(2.19) for 
the case / = 1. As discussed in |3], there are different types of solutions (generic, bouncing), 
but only one special type corresponds to confinement. The trick of obtaining this solution 
is to start numerical integration at some very large value of z or A(z) using as the initial 



condition (2.23). This automatically leads to the correct special solution. Most concretely 



this is seen from the large A limit (2.26) of the beta function; the unwanted solutions have 
/3 -^ -3A. 

Following Eqs. ( |2.23[ ) - ( |2.28[ ) we start the integration of ( |2.16[ ) - ( |2.19| ) with / = 1 using 
the initial conditions 



Hzi) -- 


= b{Xi) = 


W{zi) -- 


= WiX,) 


Kzi) -- 


= Aj, 



boXf^\UogXi 



,1/2 



l^AW) 



(A.l) 
(A.2) 
(A.3) 
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where 60 is an arbitrary constant. Now choose some numerical value for bo, Zi and Aj (say, 
bo = 1, Zi = 40, Aj = 10^'''^), and integrate the system of differential equations towards smaller 
z until the b ~ ^/{z — zjjv) divergence is encountered at the AdS boundary. Then the energy 
scale may be readjusted to a chosen value of A (as in (2.21[)) and the coordinate origin moved 



to zuv = by the scalings described in appendix A of [4J (since / = 1, only scalings 2. and 
3. are needed). This eliminates the effect of the arbitrarily chosen constants 69 and Zj. 

Note that since the approximation error is 0(1/ log A), Aj has to be very large. It also 
has to be very large in order to reach the large-z minimum Zmjn of the potential, due to the 
extremely slow increase of z(X). 

After this procedure we have a set of functions b{z), X(z) for each set of parameters c, e, 
z ranging over the interval (0,z(A = 10^'"')). Given this set, the potential Vsch(z) can be 



computed from (4.10). 

Once the potential has been computed, a number of lowest states of the eigenvalue spectruin 
can be solved using the standard shooting method. Let us choose as a reference case c = 
9/(1 + e), e = 1/100, for which 2;min.rcf = 2.06089 • 10^'^. We set in the reference case the 
boundary conditions V'(2i,r) = at z/ = 10^ and at Zr = 5 ■ 10^^, far above and below the 
minimum. For other values of c, e, we simply scale zi^r -^ -2/,r-Zmin/-2min,ref- We then solve by 
numerical rootfinding for values of m^ where the numerically computed solutions with the 
two boundary conditions can be matched such that they are continuously differentiable. 

B Numerical evaluation of the quasinormal modes 

The quasinormal modes can be numerically evaluated from the Schrodinger type equation 



(4.4); for a review, see |22j . Since the equation is linear, given a value of u; a solution is 
completely determined by two complex constants. One of them is simply a normalization, 
and the other can be determined from the in-going wave condition at the horizon, which in 
terms of u can be expressed as 

^p{u) -^ 6+*'^'' for large u. (B.l) 

This reduces the problem to finding those values of oj where the solution satisfying this 
condition is normalizable, i.e. ^(z) — >• z^ when z ~ u — t- 0. 

The standard method, used in Section [5} requires expanding the solutions around z = 0. 
However, in the present gravity dual these expansions are not power series in z but contain 



logarithmic terms, see for example (2.21). To avoid this we apply a more direct numerical 
method of removing the unnormalisable component from the solution. 

In practice, since uj is never numerically exactly at a quasinormal mode, we expect that 
the solution will contain some part of the non-normalizable (/>—)• 1 solution at the boundary, 
and V'(^) ~^ u~^''^ will therefore diverge. If we are close to a quasinormal mode, it seems 
intuitive that this divergence will happen only when u is very small (this is the heuristic of 
the algorithm, since we do not have a proof of this). This suggests a way to look for the 
quasinormal modes: given a trial lo, compute numerically the solution froin large u towards 
the boundary and find the smallest-u local minimum of |('0(w)|- Denoting its position by 
Umin, the zeros of Umin{^) are at the quasinormal modes. With finite precision numerics we 
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can then simply look for the minima of Umin{^)- See Fig. p] for a visualization of the process 
when applied to the model of Section [Sj for which the usual power expansion method could 
be applied. 



To start the integration of ( |4.4[ ) one may note Eq. (4.7), at large u 



0, 



(B.2) 



where A^'^'{zh) is as defined in Eq. (4.8). With this potential (4.4) has the analytic solution 

2 / . ,.. , . If. A „. / 2 /..,.. . i/,,A 2iuj 

V = -^. 

fh 



^{u) = Cil, 



-fh 



A^){zh)e2 



fhU 



+ C2I-, 



-fh 



A«(zft)e2^'^" 



(B.3) 



where the first term has the correct large u limit e^*"^". One can then use this evaluated 
at some u somewhat larger than Zh as the initial condition of numerical integration. This 
allows us to choose the solution with the correct asymptotic at infinity, which is the first 
term, and use the solution for the exponential potential to start the numerics at a somewhat 
smaller u. With this improvement, the algorithm works in the analytic model of Section [5] 
for Zfi = 1 ... 4, although at z/^ = 1 there still seems to be a tendency for the minima to lie 
at a too small Ima;. 



Mmin(<^) 




Figure 12: The numerical method applied to the analytic model ( |5.1[ ) with Zh = 4. Left panel: An 
example of looking for the fourth minimum. The plot shows \ip{u; lu)\ for a;„ G {4.550 — O.lOOOi, 4.573 — 
0.1083i, 4.575 - 0.1137i, 4.576- 0.1125i, 4.597- 0.1167i}- The arrows mark the location of the minima 
Umin{^)- Notice how the minima first approach the origin, then after the closest point w^ = 4.575 — 
0.1137i they again start to move away. This is the computed value of the quasinormal mode. Right 
panel: The minima Umm(w) on the complex plane. On the gray area Umin > 0.5, and Umm decreases 
for shades towards dark blue. The red dots show the positions calculated in fig. [5] for comparison. 
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